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. Abstract 
X) ■ 

. In this paper, we consider a facility location problem to find a minimum-cost coverage 

[T i ■ of n point sensors by disks centered at a fixed line. The cost of a disk with radius r has 

a form of a non-decreasing function f(r) = r a for any a > 1. The goal is to find a set of 
' disks under L p metric such that the disks are centered on the x-axis, their union covers 

the n points, and the sum of the cost of the disks is minimized. Alt et al. [T| presented 
an algorithm in 0{n A log n) time for any a > 1 under any L p metric. We present a faster 
algorithm for this problem in 0(n 2 logn) time for any a > 1 and any L p metric. 
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O ■ 1 Introduction 

■ We consider geometric facility location problems of finding k disks whose union covers a set 

P of input points with the minimum cost. A center of the disk of radius r is often modeled as 
\ a base station(server) of transmission radius r and an input point as a sensor (client), so we 

assume the cost of the disk to be r a for some real value a > 1. Thus the goal is to minimize 
r(Di) a where the disks Di covering P have radius r(-Dj). Alt et al. [T] presented a number 
of results on several problems related in this context. Among them, we focus on a restricted 
version in which the centers of the disks are restricted to be on a fixed line, simply saying 
cc-axis. When the fixed line is not given, but its orientation is fixed, finding the best line 
giving the minimum coverage even for a = 1 is quite hard to compute exactly [T], thus they 
gave a PTAS approximation algorithm. 

Alt et al. pQ presented dynamic programming algorithms for this restricted coverage prob- 
lem by aligned disks on a fixed line in time 0(n 2 logn) for a = 1, and in time 0(n 4 logn) 
for any a > 1 under any L p metric for 1 < p < oo. For metric, they presented an 
0(n 3 log n)-time algorithm. 

We reinterpret their dynamic programming algorithms together with new observations, 
then we present improved algorithms in 0(n 2 logn) time for any a > 1 and any L p metric, and 
in 0(n 2 ) time for metric. The number of disks in the optimal covering is automatically 
determined in the algorithm. If one would want to restrict the number of disks used, say as 
a fixed 1 < k < n, then we can find at most k disks whose union covers the input points 
with minimum cost in a similar way. Actually we can find such k disks for all 1 < k < n in 
0(n 3 logn) time in total. 

The formal definition of the problem is as follows: Given a set P = {pi,P2, ■ ■ ■ ,Pn) of n 
points in the plane, a real value a > 1 and L p metric for some p > 1, find an optimal disks 
Di, Z?2, • • • , Dk with centers Si on the x-axis and with radii r{Di) whose union covers P such 
that the sum of the radii, ^2ir a (Di) is minimized. 
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2 Geometric properties 



We assume that the line where the centers of the disks lie is x-axis. As mentioned in [?], 
we assume that all points in P lie above or on the x-axis and no two points have the same 
x-coordinates. If a point p is below the x-axis, we replace it with a new point p' mirroring 
p with respect to the x-axis, then we get the same optimal covering. If p is directly above 
p' , then any disk containing p always contains p', so we can simply discard p from P. Thus 
from now on we assume that the points of P have nonnegative and distinct x-coordinates, 
and they are indexed from left to right. Finally we assume the points of P are in the general 
position, i.e., no three or more points lie on the boundary of a disk with centers on the x-axis. 

We also notice that the optimal covering is not unique, so we assign the lexicographic 
order to the optimal covering, the set of the disks according to x-coordinates of their centers. 
Then we consider only the leftmost optimal covering D = {Di, Di, ■ ■ ■ , D^} with centers in 
increasing order on x-axis. 

Let a < 1 and let r(Z)j) denote the radius of Di. We call r a {Di) the cost of the disk Di. 
For a while, let us consider L p metric only for 1 < p < 00. Let dR denote the boundary of 
a closed region R. We denote by U the highest point(or apex) of dDi, and by ai and hi the 
left and right intersection points of dDi with the x-axis, respectively. Let B be the union of 
disks in T>. Then the following facts hold; the first one is mentioned also in pQ. 

Fact 1 [1] For each 1 < i < k, the apex ti of Di appears on dB. 

Let us consider dDi n dB, i.e., the circular arc of dDi which appears on dB. By Fact [H ti 
must be contained on the arc, so the arc is divided into the left and right subarcs at ti. Then 
we have the following fact. 

Fact 2 For each 1 < % < k, dDi n dB must contain either one point of P at the apex ti or 
two points of P, one on the left subarc and the other on the right subarc of dDi CI dB. 

Proof. It is obvious that there must be at least one point of P on dDi n dB. Otherwise 
we can shrink Di to get a smaller cost until dDi contains some point. Also if one of the left 
and right arc has no points, then we can shrink Di while keeping the point on the one subarc 
until some point lie either on the apex ti or on the other subarc containing no points. This 
contradicts to the optimality. EO 

For each 1 < i < k, we define £i as a vertical line between Di and -D^+i; if Di intersects 
Di + i, then li is a vertical line through intersections dDi n dDi + \, otherwise li is an arbitrary 
vertical line between bi and dj+i. For convenience, we define £0 and as vertical lines passing 
through a\ and bk, respectively. 

Let Pi be a subset of points of P lying between and li for 1 < i < k. Then we know 
that Pi contains at least one point by Fact [2j and they are pairwise disjoint and their union is 
the same as the whole set P. Let Cj be the smallest axis-centered disk containing P,. Clearly 
{Ci, . . . , Ct\ is a covering for P. We have the following lemma. 

Lemma 1 Ei^^Q) = £i<i< fc r«(A). 

Proof. Since {D\, . . . , D^} is the optimal covering for P, it holds that £V r a {Di) < ^ r a (C{). 
For each 1 < i < k, Pi is contained in P D Di, thus r a (Ci) < r a {Di). Since /(r) = r a is a 
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Figure 1: The farthest Voronoi diagram for P partitions the x-axis into intervals. 



nondecreasing function for a > 1, ^2 i r a {Ci) < ^2if a (Di), which completes the lemma. 1=0 

The above lemma means that there is a vertical partition of P into Pj's such that the 
smallest disks containing Pj's are the optimal disks for P. Using this lemma, we can derive a 
fast dynamic programming algorithm. 

3 Dynamic programming algorithm 

Alt et al. p] defined a pinned disk(or circle) as the leftmost smallest axis-centered disks 
enclosing some fixed subset of points, so the pinned disk contains at least one point on its 
boundary. The disk Cj defined in Lemma Q] is a pinned disk. It is obvious that the optimal 
covering P is a subset of such pinned disks. In [I], the dynamic programming algorithm 
chooses pinned disks with minimum cost from all 0(n 2 ) pre-computed pinned disks, satisfying 
the feasibility condition that no other points of P lie above the chosen pinned disks. This 
step causes to take the total 0(n 4 log n) time. But Lemma [1] tells us there must be a partition 
Pi, . . . , Pk, separated by vertical lines, such that a set of the smallest disks containing Pi is 
indeed an optimal covering for P. Thus we simply go through the input points from left to 
right, not through the pinned disks, and compute the smallest disk Cj enclosing Pj instead of 
checking the feasibility condition. 

Let A be an array in which A[i] stores the minimum cost for a subset {pi,Pi+i, . . . ,p n }. 
The minimum cost for the whole set {pi, . . . ,p n } will be stored at A[l]. If we denote by 
D({pi, . . . ,Pj}) the smallest disk containing {pi, . . . ,pj}, then we have the following recurrence 
relation: 



The key step is to compute D({pi, . . . ,Pj}) fast. We can do this in amortized O(logn) 
time maintaining the intersection of the x-axis with the farthest Voronoi diagram(FVD) in a 
dynamic way. For a fixed i, A[i] is computed in O(nlogn) time, so the total time to compute 
A[l] becomes 0(n 2 logn). 




3 



As in Figure [H the intersection of the farthest Voronoi diagram for {pi , . . . ,pA with the 
x-axis partitions the x-axis into intervals I\, I2, ■ ■ ■ , I m (i,j) from the left to the right, where Ii 
is a half-open interval // := [x;_i,x;), where x = —00 and x m Uj\ = +00. Each interval // is 
a collection of the points from which the farthest point of {pi, . . . ,pj} is the same. We denote 
by p(I\) the farthest point from any x € 1\. Then a disk centered at some point x € and 
with radius \xp(I{)\ encloses all the points of {pi, . . . ,Pj}- 

Let D(Ii) be the smallest disk enclosing {pi, . . . ,pj} whose center lies in I\. We have two 
cases. For a case which dD(I[) has one point at its apex, the point is indeed p{Ii) and the 
center of D{I{) has the same x-coordinate as that of p{h)- For the other case, dD(Ii) should 
have two points, so the center of D(Ii) must be on x/_i, the left endpoint of I[, but the radius 
of D(Ii) is defined as 00. 

To store such intervals, we use a balanced search tree T [3]. We store at its leaves the 
intervals I±, . . . , I m aj) with their corresponding radii from left to right. Each internal node v 
of T stores the minimum one among the radii in the leaves of the subtree rooted at v. Then 
the radius stored at the root of T is the radius of the smallest disk enclosing {pt, . . . ,Pj}. We 
can insert a new interval into T and delete an interval from T both in O(logra) time. 

For a fixed i, we now construct the intervals I±, . . . , I m u } j) for all i < j < n incrementally 
from j = i to j = n. For j = i, there is only one interval. We start with this interval, and 
update the interval set by adding the points one by one from pi+i, . . . ,p n . We now explain 
how we update T for {p,, . . . ,Pj~i} when pj is inserted. 

We know that the interval for pj must appear because pj is the rightmost point among 
Pi, . . . ,Pj-i, and moreover the interval should be the leftmost one, i.e., its left endpoint must 
be xo = —00. When the interval for pj is inserted into T, several consecutive intervals in T 
from the left should be removed from T or replaced with a shorter interval in T. To identify 
such intervals, we need the following basic properties on the farthest Voronoi diagram. 

Lemma 2 For {pi, . . . ,pj} under any L p metric, 1 < p < 00, the intersection of the x-axis 
with the farthest Voronoi diagram for {pi, . . . ,pj} has the properties: (1) The interval for pj 
is connected, and (2) for any two consecutive intervals I and J where I is in the left of J on 
the x-axis, then p{I) > p(J), where p(I) > p{J) means the x-coordinate of p(I) is larger than 
that of p{J). 

Proof. For the completeness, we prove these properties. A bisector of two points under any 
L p metric is monotone to the x-axis and the y-axis, so it intersects the x-axis only once [4]. 
To prove the connectedness, we suppose that pj has two disjoint intervals / and L, where / is 
to the left of L. There must be one or more intervals between them, denote by J the interval 
to the right of / and by K the interval to the left of L. Note that J is not necessarily different 
with K. Let D be a smallest disk centered at I n J, i.e., the common endpoint of I and J 
which encloses all points in {pi, . . . ,Pj}- Then p(I) and p{J) lie on dD. Similarly, let D' be a 
smallest disk centered at KnL enclosing all the points. Since p(I) = p(L) = pj, they must be 
on one of two intersections dD n dD', clearly the one above the x-axis. Also the lune D n D' 
contains all the points in {pi, . . . ,Pj}. This implies that p{J) must lie on the right boundary 
arc of the lune. The bisector of p(I) and p{ J) intersects the x-axis at I C] J, thus the points 
on the x-axis to the left of / n J is farther to p(J) than to p(I), which contradicts that I is 
in the left of J. For the second fact, we consider the half-circle of the smallest disk centered 
at In J on the x-axis which passes through p(I) and p{J). Since the half-circle intersects 
with the bisector of p(I) and p(J) exactly once, p(I) should be in the right of p(J) along the 
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half-circle. This means p(I) > p{J) because the half-circle is monotone to the x-axis. 

Let J = [a, b) be the interval of pj in the interval set for {p^, . . . ,pj}. Then we already 
know that a = -co. By Lemma [21 it suffices to find the interval from the intervals for 
{pi, . . . ,pj_i} which intersects with the bisector of pj and p{Ii). Then b is the intersection 
of Ii with the bisector. For this, we do the intersection test from I = j — 1 to I = i one by 
one. Once Ii is found, we (1) delete the intervals Ii, ■ ■ ■ , Ii-i, which are completely contained 
in J, from T, (2) insert a new interval J for pj, and (3) replace(i.e., delete then insert) I\ 
with a part not contained in J, l\ \ J. If some interval is removed from T, then it is never 
inserted again into T. Hence, for a fixed i, we can compute the smallest disks enclosing disks 
for {pi, . . . ,pj} for all i < j < n in 0{{n — i) log n) time. In other words, we can compute A[i] 
in 0((n — i)logn) = 0{n log n) for fixed i. The total time of the algorithm is 0(n 2 log n), 
and the space is 0(n). The detailed algorithm is summarized below. 

Algorithm 1 MinCostAlignedCoverage(P, a) 
Input: A set P of n points {pi, . . . ,p n } and a > 1. 

Output: A set of disks T> = {D\, . . . , D^} with minimum cost of r a {Di) which covers P. 
l: A[n + 1] = oo. 

2: Initially, T consists of one interval (— oo, +oo) for p n 

3: for i n to 1 do 

4: A[i] = OO 

5: for j <— i to n do 

6: Find the first interval // in T such that the bisector B of pj and 

intersects 1\ by scanning the intervals in T one by one from left to right 
7: J := [-oo,5n It) 

8: Remove intervals Ii, . . . , Ii-i, replace /; with I\ \ J, and insert J in T 
9: Let r be the radius stored at the roof of T, i.e., r = r(D({pi, . . . ,Pj})) 

10: A[i] = mm(A[i], r a + A[j + 1]) 

11: Keep the index j which gives the minimum cost 

12: end for j 

13: end for i 

14: Reconstruct the optimal disk set T> by backtracking the recorded indices 
15: return A[l] and T> 



Algorithm for metric. Under this metric, the unit disk is an axis-aligned square. As 
before, we consider only the leftmost optimal covering by the lexicographic order. We can 
easily see that Fact Q] and Fact [2] can be applied for L^, metric if ti, the apex of the disk is 
defined as the upper and right corner of the disk. To use Lemma O we define a partition of 
P, Pi,... ,Pk, separated by vertical lines containing right sides of the optimal disks. Then 
we can also prove in a similar way as the proof in Lemma [2] that the sum of the costs of the 
smallest squares C{ containing Pi is the same as the minimum cost for P. We now compute 
A[i] similarly. The key step is to compute the smallest square C enclosing {pi, . . . ,pj} quickly. 
This square C is determined by two points; pj and one of the points pt and the highest point 
of {pi, . . . ,pj-i}, which can be computed in O(l) time if we maintain the highest point during 
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the incremental evaluation. Thus we can compute A[i] in 0(n) time. The total time is 0(n 2 ). 

Theorem 1 Given a set P of n points in the plane and a non- decreasing cost function with 
a > 1, we can compute an optimal disks centered on the x-axis such that the union covers P 
and the sum of the costs of the disks is minimized in 0(n 2 log n) time for any fixed L p metric 
and in 0(n 2 ) time for metric. 

We can also consider the case when the number of disks used to cover P is given as a 
fixed value k. This case would be required by practical reasons. This can be similarly solved 
by filling a two dimensional table ^4[i][/c], the minimum cost needed to cover Pi,...p n with 
at most k disks, in 0(kn 2 logra) time. Actually we can find all optimal coverings for any 
1 < k < n in the same time. 

Theorem 2 Given a set P of n points in the plane and a non- decreasing cost function with 
a > 1, we can compute a collection of all optimal coverings for P such that P is covered by 
at most k disks for any 1 < k < n and the sum of the costs of the disks is minimized in 
0(n 3 log n) time for any fixed L p metric and in 0(n 3 ) time for L ro metric. 

4 Concluding Remarks 

We can consider other disk coverage problems with practical restrictions such as the connec- 
tivity constraint. Recently, Chambers et al. [2] investigated a problem of assigning radii to a 
given set of points in the plane such that the resulting set of disks is connected and the sum 
of radii, i.e., a = 1 is minimized. When we bring such connectivity constraint to our problem 
for a > 1, we need to find a "connected" set of disks centered on the x-axis whose union 
covers n input points. When a = 1, the smallest disk containing all points is the optimal 
coverage. However, we can easily show for a > 1 that infinitely many disks always guarantee 
the minimum cost coverage for any input. Thus we should restrict the number of disks used 
to cover, say 1 < k < n. But we have no idea how hard this problem is for a fixed k. 
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